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SECTION I. INTRODUCTION 


This report sununarizes the work vinder Contract N?VS8-33433 
during the 12-month period beginning April 10, 1980. The de- 
tailed background information concerning earlier work on this 
contract can be found in References 1 and 2. The objective of 
the study is the development of a vector wind gust model that 
is suitable for orbital flight test operations and trade 
studies. During this reporting period, emphasis has been given 
to verification of the hypothesis that gust component variables 
are gamma distributed, gust modulus is approximately Weibull 
distributed, and zonal and meridional gust components are 
bivariate gamma distributed. The body of the report is con- 
tained in four sections (II through V) . Section II describes 
a method of testing for bivariate gamma distributed variables; 
in Section III, two distributions for gust modulus are described 
cind the results of extensive hypothesis testing of one of the 
distributions are presented; Section IV establishes the validity 
of the gamma distribution for representation of gust component 
variables. Conclusions are presented in Section V. 


SECTION II. TESTING FOR BIVARIATE GAMMA 
DISTRIBUTED VARIABLES 


The hypothesis that absolute component gust and associ- 
ated gust length are bivariate gamma distributed can be tested 
according to the procedure described below. 

Consider variables x and y that are distributed ac- 
cording to the bivariate gamma distribution with scale param- 
eter 0 ^ and @2 • dimensionless variables and T2 are 

defined by 



( 1 ) 


The variables Tj^ and T2 can be expressed in a coordinate 

system that is rotated by 45 *; the transformed variables 
Zj^ and 22 are given by 
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It can be shown that the probability, P^, that bivari- 
ate gamma distributed variables and Z2 will occur within 

the area bounded by the lines z^ - Z2» and z^^ = -Z2 

(illustrated in Figure 1 ) can be estimated from the series: 
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Figure 1. Area, A (Shaded) , Which Bounds Bivariate Gamma 
Distributed Variables and Z2 for Which a 

Probability of Occurrence Can Be Calculated 
from Equation 3 
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H(a,X) is the incon^lete gamma function which is given 
by the series 


H(a,X) 


a “X " X^ 

* „!o TTs+imr 


( 4 ) 


where 


a = 2(y+m) 


p = correlation between variables x and y 


Alternatively, can be calculated by numerical 


integration of the equation 


m fH 


-/2z. 




= 




(1-p)^ r(Y) 


(5) 


where I (r^ z, } is the modified Bessel function of the first 
n l-p i 

kind of order n; for Y * Zy^, I^(Y) is calculated with the 
series approximation, 

„n+2k 

I_(Y) = I 

" k=0 2*^ ‘^*^kJ r(n+k+l) 


A computer program for calculation of P^, using 

Equation (3), has been developed. A sample of the calculations 
of P^ as a function of p for y= 3 are listed in Table 1. 
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Table 1. Pa (P » Y=3) Calculated According to Equation 
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standard results for checking the computer programs are 
obtained from the closed form solution for the case y“l which 
can be expressed in the form 

Y 

7p 

P^(Y=1,P) - 1 - [e^{^ + 1) - e“^(^ - 1)1 (7) 


Values for are listed in Table 2 for selected values 

of zj and p. 


Table 2. P^(y» 1, p) Calculated from Equation (7) 


1 . 

P 

N 



.1 

.25 

.50 

.75 

.875 

1 

.425980 

.445255 

.474470 

.495090 

.501805 

2 

.774586 

.774144 

.769772 

.763367 

.760084 

3 

.919308 

.911445 

.899446 

.889099 

.884464 

4 

.971975 

.965471 

.956084 

.948025 

.944361 

5 

.990370 

.986548 

.980820 

.975641 

.973206 

6 

.996704 

.994760 

.991624 

.988584 

.987097 

7 

.998874 

.997959 

.996342 

.994650 

.993786 

8 

.999615 

.9®9205 

.998402 

.997493 

.997008 

9 

.999869 

.999690 

.999302 

.998825 

.998559 

10 

.999955 

.999879 

.999695 

.999449 

.999306 


Another useful special case is for p=0 and 2 y equal to an 
integer. 

P^(p=0, 2 y * an integer) 


= 1 


/2 zj 


2y-1 

I 

k=0 


2k/2jz*)’' 

in 


( 8 ) 
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The variation of as a function of correlati^a 

coefficient, p, (for y*2) and as a function of shape parameter, 

Y, (for p>.5) is illustrated in Figures 2 and 3, reopectively . 

A comparison of observed and expected is illustrated 

in Figure 4; the line drawn at an angle of 45* to the abscissa 
represents perfect agreement between observed and expected 
values. Deviations of the plotted points from the line repre- 
sent differences between the observed and expected values. The 
data plotted in Figure 4 show a consistent pattern at 10 and 
12 km; for P^ < .3, the observed is larger than the expected; 

for intermediate values (.3 < < .8), the expected is larger 

than the observed. These results are the basis for initiating 
a more detailed analysis of the validity of the gamma distribution 
hypothesis for the marginal distributions (component gust and 
associated gust length) . The results of this analysis are 
described in the next section. 
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Figure 3. Series Approximation of P* as a Function of z* and y for 






























SECTION III. THE OISTRIBUTION OF GUST MODULUS 


Given that the absolute gust components are uncorrelated 
bivariate gamma distributed, it follows that the probability 
distribution of the gust modulus, R, can be calculated with 
the double series approximation. 


PrCH 1 R*} = G(R*) = 


Y1+Y2 ^ 

R* ^ 02 

~~irVf[TTTy^ 


z 

i=0 


(-R* 02 >^ 


/Y1+Y2+M 


i 

• Z 

n =0 


0 ^" 02 "”r(~ 


Yi+n, 


(V)r( 

n ! (i-n) ! 


Y2+i-n 


( 9 ) 


where and 3 2 scale parameters and y and Y2 the 

shape parameters of the u and v component gust distributions, 
respectively. 

Smith* has shown that the above expression is approximately 
equivalent to 


G(R*) 


H(Yj^ + Y2» AR*) 
+ Y2) 


( 10 ) 


where H(Yj^ + Y2» AR*) is the incomplete gamma function which can 

be calculated accurately with the series approximation given in 
Section II. 


*Personal communication. 
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A = 


( 11 ) 


Yi Y, ^1 ^2 

r{-f) T(-f) r(Yi+Y2) ®2 

/Y 1 ■♦’Y 2 \ 

2r(Y^) r(Y2)r(-V^j 


1 


Preliminary tests have indicated that reasonably accurate 
estimates of the probability distribution can be obtained from 
equation 10. However, it would be advantageous to determine 
if an alternative expression can be found which would not re- 
quire as much computation. The Weibull distribution, widely 
used in wind energy studies (Reierence 3) was chosen to repre- 
sent gust modulus because of its relative mathematical simplicity 
and the availability of data for parameter estimation. The cumu- 
lative probability function for the Weibull distribution of 
gust modulus is 


G(R*) = 1 - EXP 


c 


( 12 ) 


The parameters k and c are calculated according to the 
approximation given by Justus (Ref. 3) 


/ Oj^\ -1.086 

\ R/ 


(13) 


rTT + rm 


(14) 


It is noted that equation 13 implies the relation. 


/IbV = j^-1. 84162 
\ R / 


(15) 
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whereas the exact relation for a Weibull distribution is 
given by 



r(l -I- 2/k) 

/ 

J^d + 1/k) . 


“ 1 


( 16 ) 


The accuracy of the approximation has been evaluated for 
values of k from 0.5 to 10 by calculating the ratio, P, of 
the right side of equation (16) to the right side of equation (15) 
Perfect agreement is indicated when P=l. As illustrated in 
Figure 5, for k > 1, P is within a few percent of unity; for 
k < 1, P approacFes » as k approaches 0. Therefore, it is con- 
cluded that the approximation given by equation (13) is accurate 
for k ^ 1. The calculated values of k for gust modulus are 
between 2 and 3 which is within the range of acceptable accuracy 
of equation (13) . 

A comparison of the Weibull, The probability distribu- 
tion associated with the modulus of a bivariate normal, and 
the observed probability distribution is illustrated in 
Figure 6. It is indicated that there is little difference 
between the theoretical distributions for percentiles between 
20 and 98; for percentiles outside that range, the distribu- 
tions diverge; for this case, the observed distribution fits 
the Weibull slightly better than the bivariate gamma modulus 
distribution. 

The hypothesis that gust modulus at a reference altitude 
is drawn from a Weibull distributed population was tested for 
69 cases. The results summarized in Table 3 indicate that the 
hypothesis is accepted at the .05 level of significance in a 
large majority (65/69) of the cases. 
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Figure 6. Observed and Theoretical Distribution of Gust Modulus at 12 km 
During February at Cape Kennedy for X = 2,470 m 



Table 3. Summairy of Results of Testing the Hypothesis* That 
Gust Modulus at a Reference Altitude (4, 6 ... 14 km) 

Is Drawn From a Weibull Distributed Population 
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SECTION IV. DISTRIBUTION OF GUST COMPONENT VARIABLES 


Four variables associated with gusts at a reference height « 
Hq, have been studied to establish the validity of the hypothesis 

that they are samples from gamma distributed populations. The 
four variables are illustrated in Figure 7. The variable u^ is 

the largest u component excursion with sign equal to the sign of 
u at Hq; U 2 is the largest u con^onent excursion of sign opposite 

Uj^ found by scanning upward after the second zero crossing associ- 
ated with Uj^. The vertical distance between and U 2 is defined 
as L Range; the sum of the absolute values of and U 2 is defined 

as u Range. The variables u Range and L Range represent wind 
shear and wind shear altitude interval associated with gusts in 
the vicinity of Hg. Each of the four variables defined above have 

been calculated at six reference altitudes from a sample (150/ 
month) of February, April, and July Jimsphere wind profile data 
from Cape Kennedy. These data sets were tested to establish the 
validity of the hypothesis that each variable is drawn from a 
gamma distributed population. Acceptance or rejection of the 
hypothesis is at the .05 level of significance for a variate 
defined by 


X 


2 


N (O^ 
i=l 



2 


(17) 


= Observed frequency in the ith class interval 

E. = Expected frequency in the ith class interval 
(of the theoretical -gamma distribution) 


The results of the hypothesis testing are described below. 


A, Absolute Gust Component and Associated Gust Length 


Gust, defined as the maximum excursion between successive 
zero crossings in the vicinity of a reference altitude, and 
associated gust length, defined as the distance between zero 
crossings, are each hypothesized to be drawn from a gamma 
distributed population. The hypothesis is accepted at the 
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.05 level of significance^ in a large majority of cases/ for 
gust component magnitude (|u'|); specifically/ the accept/ 
reject ratio is 47/22 and 46/23 for u and v coR^onent magni- 
tudes/ respectively. As indicated in Table 4/ the ratio is 
significantly smaller for gust length (Lu and Lv) with rejec- 
tions exceeding acceptances (for method I). The large number of 
rejections is attributed to large differences between observed 
and expected frequency of occurrence in the first few class 
intervals; the observed frequencies are always much larger than 
the expected frequencies. Small gust magnitudes are associ- 
ated with small gust lengths that are observed as a consequence 
of the definition of gust used in this study. These small gust 
lengths are not measurable with the Jimsphere system; therefore/ 
they are not considered to be valid data for hypothesis testing. 
By neglecting these data, we obtain the results summarized 
under II in Table 4 which indicate acceptance in a much larger 
proportion of the cases. 


B. U Range and L Range 


A summary of results of testing the hypothesis that the 
variables, U range and L rcmge, are drawn from gamma distributed 
populations is given in Table 5. It is indicated that the 
hypothesis for U range is accepted at the .05 level of signifi- 
cance in 66 of the 72 cases. Acceptance is not a function of 
altitude except in July when the number of samples accepted at 
14 km was less than at the other altitudes. Acceptance was 
not related to filter choice with only slight exceptions 

(for A * 2470 during July and A = 6000 m during February 

c c 

one-third of the samples were rejected). Based on these results/ 
it is concluded that U range is gamma distributed. 

The results for L range summarized in the lower half of 
Table 5 indicate acceptance of the hypothesis (46 of the 72 
cases) with not as strong a tendency as that indicated above 
for U range. Acceptance is an irregular function of altitude 
which is a minimum at 12 km where 50 percent are accepted to 
a maximum at 8 km where 75 percent are accepted. Acceptance 
is greater in July (75 percent) than in either April or 
February (58 percent for both months) . Acceptance is weak or 

non-existent for A * 420 m and is strong for A large (2470 

c c 

and 6000 m) . 
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Table 4. Summary of Results of Testing the Hypothesis' 
that u and v Component Absolute Gust and Gust Length are 
Drawn from Gamma Distributed Populations 



All Filters 


Filter 

A (m) 



420 

997 

2470 

6000 

6/0 

4/2 

3/3 

3/2 

All Filters 

16/7 



All 

Months 

2 

16/2 

6/0 

11/7 

5/1 

10/8 

4/2 

9/6 

5/0 

46/23 

20/3 


All 

Months 




420 

5/1 

4/2 

997 

0/6 

2/4 

2470 

1/5 

3/3 

6000 

2/3 

1/4 

All Filters 1 

8/15 

10/13 



2/6 

5/1 

5/1 

5/1 

15/3 

4/14 

1/5 

5/1 

4/2 

10/e 

7/11 

3/3 

4/2 

5/1 

12/6 

7/8 

4/1 

4/1 

'■ ■ ■ 1 

4/1 

12/3 



^^^At the .0$ level of significance for variate with » degrees of freedoat n • n*l-b, tdiere 
n • rusaber of class intervals « b • n%mbcr of paraaeters of g sMS i distribution ■ 2. 

iS the ratio of the nunber of cases accepted to the ntasber rejected. 
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Table 5. Summary of Results of Testing the Hypothesis 
That the Variables, U Range and L Range, at a 
Reference Altitude (4, 6, ... 14 km) are 
Dravm from Gamma Distributed Populations 


mt«r 


Raferance 

Altltuda 

On) 


Inssary 

Veriah: Month 

4 

6 

i 

10 

12 

14 

A R 

U range 2 420 

A* 

A 

A 

A 

A 

A 

6 0 

997 

A 


A 

A 

A 

A 

6 0 

2470 

A 

h 

A 

A 

A 

A 

6 0 

6000 

R* 

A 

k. 

A 

R 

A 

4 2 

Accept/Reject (all filters) 

3/1 

4/0 

4/0 

4/0 

3/1 

4/0 

22 / 2 

4 420 

A 

A 

A 

A 

A 

A 

6 0 

997 

A 

A 

A 

A 

A 

A 

6 0 

2470 

A 

A 

A 

A 

A 

A 

6 0 

6000 

A 

A 

A 

A 

R 

A 

5 1 

Xccept/Ra ject 

4/0 

4/0 

4/0 

4/0 

3/1 

4/0 

23 / 1 

7 420 

A 

A 

A 

A 

A 

A 

6 0 

997 

A 

A 

A 

A 

A 

R 

5 1 

2470 

A 

A 

A 

R 

A 

R 

4 2 

6000 

A 

A 

A 

A 

A 

A 

6 0 

Accept/Raject 

4/0 

4/0 

4/0 

3/1 

4/0 

2/2 

21 / 3 

Accept/Reject (all months) 

11/1 

12/0 

12/0 

ll/l 

10/2 

10/2 

66/6 

L range 2 420 

R 

R 

R 

R 

R 

R 

0 6 

997 

A 

A 

R 

A 

R 

R 

3 3 

2470 

A 

A 

A 

A 

A 

A 

6 0 

6000 

R 

A 

A 

A 

A 

A 

5 1 

Acccpt/Re ject 

2/2 

3/1 

2/2 

3/1 

2/2 

2/2 

14 /lO 

4 420 

R 

R 

A 

R 

R 

R 

1 5 

997 

A 

A 

R 

A 

R 

R 

3 3 

2470 

A 

A 

R 

R 

A 

A 

4 2 

6000 

A 

A 

A 

A 

A 

A 

6 0 

Accept/Pe ject 

3/1 

3/1 

2/2 

2/2 

2/2 

2/2 

14 /lO 

7 420 

R 

R 

A 

R 

R 

R 

1 S 

997 

A 

A 

A 

A 

R 

A 

6 1 

2470 

A 

A 

A 

\ 

A 

A 

6 0 

6000 

A 

A 

A 

A 

A 

A 

6 0 

Accept/Ro ject 

3/1 

3/1 

4/0 

3/1 

2/2 

3/1 

16/6 

Accept/He jcct (all months) 

3/4 

9/3 

8/4 

8/4 

6/6 

7/5 

46/ 26 

•Accept (A) or Reject (K) hypothesis 

at the . 

,0S level 

of significance 

for 

variate 


with m degrees of freedom, m • n-l-b, where n • number of dess intervals, b • miabcr 
of {idr^eters of the ganrui distribution “ 2. 


SECTION V. CONCLUSIONS 


This report has en^hasized methods for establishing the 
validity of the hypothesis that observed gust variables, 
iiiCluding gust component magnitude, gust length, u Range, and 
L Range, have been drawn from gamma distributed populations 
and that observed gust modulus has been drawn from a bivariate 
gamma distributed population that can be approximated with a 
Weibull distribution. An analytical procedure has been pro- 
posed for testing for the bivariate gamma distribution. The 
procedure has the advantage of not requiring frequency counts 
within narrow cells defined by the intersection of intervals 
of the marginal distribution; these frequency counts would be 
impractical and unreliable because of the limited sample size 
(150) of the available data. Instead, the new method requires 
theoretical and observed frequency counting over larger areas 
associated with non-dimens ionali zed and transformed variables. 
Preliminary results utilizing this method have indicated larger 
observed than expected frequencies for small gust lengths and 
associated small gust magnitudes; this is attributed^le to the 
definition of gust used in this study. These small gust lengths 
are not measurable with the Jimsphere system and as indicated 
in Section IV the results of hypothesis testing for the margi- 
nal distributions are improved greatly by eliminating them from 
the data sample. The hypothesis that gust con^onent (u and v) 
magnitudes are drawn from a gamma distributed population is 
accepted at the .05 level of significance in 122 of the 136 
cases tested; for gust length (Lu and Lv) , 102 of the 136 cases 
were accepted. 

The variables u Range and L Range have been used to repre- 
sent component wind shear and shear interval associated with 
gusts. The hypothesis that u Range observations were drawn 
from a gamma distributed population was accepted at the .05 
level in 66 of the 72 cases tested; the acceptance ratio was 
somewhat smaller for L P^nge with acceptance in 46 of the 72 
cases tested. 

Testing of the hypothesis that gust modulus is drawn from 
a Weibull distributed population has yielded highly favorable 
results with acceptance of the hypothesis at the .05 level in 
65 of the 69 cases tested. 
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APPENDIX 


Parameters y and 6 (calculated from sample moments) for 
hypothetical gamma distributions of gust component variables 
|u'|, Iv'l, Lu, Lv, u Range, and L Range defined in Section IV 
are listed in Tables A-1 through A-6. 

The parameters in the tables can be used to derive 
the gamma probability density function of the form 


g(x) 


rTTr 


Y-i 


EXP(-Bx) 


( 1 ) 


Equation (1) 
variable y , i.e. , 


can be expressed in terms of a nondimensional 
y = ^ , such that 

P 


g(y) = EXP(-y) 


(2) 


The probability that y does not exceed a specified value, 
Y, is given by 


y 

P {y ^ y) = / g(y) dy 
0 


TOT 



0 


EXP(-y)dy 


(3) 


The integral on the right side of Equation 3 is the incom- 
plete gamma function, H(y,Y), which can be approximated with 
the series summation given by Equation 4 in Section II with 
the substitution 

a = Y 


A-1 


X = Y 


Table A-.l. Gamma Distribution Parameters y and 6 of Absolute u Component Gust 
Estimated from Sample Moment Statistics* 
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Table A-2. Ganna Distribution Parameters y 0 of Gust Length, Lu 

Estimated from Sample Moment Statistics* 
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Table A-4. Gamma Distribution Parameters y an<3 6 of Gust Length, Lv 
Estimated from Saitple Moment Statistics* 
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Table A-5. Gamma Distribution Parameters y and 6 of u Range Estimated from Sample Moment Statistics* 
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Table A-6. Gamma Distribution Parameters y and 

























